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Baryon Acoustic Oscillations (BAO) have recently been observed in the distribution of distant 
galaxies. The height and location of the BAO peak are strong discriminators of cosmological pa- 
rameters. Here we consider the ways in which weak gravitational lensing distorts the BAO signal. 
We find two effects that can affect the height of the BAO peak in the correlation function at the 
percent level but that do not significantly impact the position of the peak and the measurement of 
the sound horizon. BAO turn out to be robust cosmological standard rulers. 

PACS numbers: 98.80.Es, 98.62.Py, 98.65.-r 



I. INTRODUCTION 

In the early universe, before the recombination of pro- 
tons and electrons into neutral hydrogen, the photons, 
electrons, and protons were tightly coupled and behaved 
like a single fluid. The density of this fluid underwent 
acoustic oscillations until recombination. These oscilla- 
tions left their imprints in the spectrum of anisotropics 
in the cosmic microwave background (CMB) as well as 
in the distribution of matter via the correlation func- 
tion of the galaxy distribution. The effect in the former 
case is quite pronounced because the photons have es- 
sentially traveled freely since their last scattering at re- 
combination. In the latter case, called baryon acoustic 
oscillations (BAO), the effect is diluted since the dom- 
inant dark matter did not participate in the primordial 
dance P, H, d, 13] ■ Observations of both of these effects 
over the last few years 0, Q provide dramatic proof that 
our picture of the early universe is consistent. 

We are now free to go further and use these observa- 
tions to constrain cosmological parameters. In the case of 
BAO, the correlation function should peak at a comoving 
scale of order rp ~ 100/i~^Mpc where the Hubble con- 
stant is parameterized as Hq = lOOh km sec~^ Mpc~^. In 
a survey of galaxies at the same redshift, the peak shows 
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up at an angular scale S9 = rp/[{l + z)Da{z)] where 
Da{z) is the angular diameter distance out to redshift 
z. In a ID radial survey with redshifts, the peak will 
show up at Sz — rpH{z). More generally, in a 3D survey, 
one can hope to measure both Da{z) and H{z). These 
quantities are enormously important to cosmologists be- 
cause they depend on the energy density back to redshift 
z and therefore on the properties of matter and dark en- 
ergy [3, H, H, ■ Besides its location, the height of the 
BAO peak is governed by the matter density VLmh? pTl |. 

For these reasons, a number of ambitious future 
surveys have been proposed aiming to measure the 
BAO peak at multiple redshifts with very high ac- 
curacy d, Roughly, one expects to measure 
Da from the acoustic scale with an accuracy Q 
{V/?>h-^ Gpc)-i/2(fc„iax/0.2/iMpc-i)-i/2 where V is the 
volume of the survey, and fcmax is the comoving wavenum- 
ber at which the power spectrum peaks. Currently for 
the Sloan Digital Sky Survey (SDSS), the error is of or- 
der 4% at z ~ 0.4 [O]. These errors are expected to go 
down [11] to 2.8% at z - 0.4 by 2008 when SDSS has 
accumulated 8000 square degrees; tol.l%atz~0.6 and 
1.4% at z - 2.5 (SDSS-III, 10000 deg^); 1% at z ~ 1 and 
1.5% at z - 2.8 (SUBARU+NOAO, 2013); and 1.2% 
at z - 0.4,0.6,1 (GWFMOS, 2000 deg^). Amid this 
excitement, many studies have been carried out search- 
ing for sources of systematic errors in BAO measure- 
ment [l3. [isl. [l6l| focusing mostly on the non- linear evolu- 
tion of structure, scale-dependent bias and errors in the 
survey window function estimation, the first two leading 
to effects smaller than a percent on Da while the latter 
has to be kept below 2% not to bias the acoustic scale by 
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more than 1% at 2: = 1. Given the great wealth of data 
that will become available, the present analysis focuses 
mainly on galaxy surveys and on the resulting galaxy 
correlation function. It is necessary to stress, however, 
that the same kind of analysis can be applied equally well 
to the two point correlation function of different classes 
of objects, notably the 10^ high redshift QSO that are 
expected to be surveyed by SDSS-III. 

In this paper we show that weak lensing is not nec- 
essarily negligible at the level of accuracy that will be 
attained by future surveys. This effect is reminiscent of 
the weak lensing of Type la supernovae that induces an 
external dispersion in their observed brightnesses, com- 
parable in magnitude to the intrinsic dispersion for red- 
shifts z > 1 [17, 18, 19, 20]. In that case, the main effect 
of lensing is to magnify the background objects; in this 
case, while magnification can cause one effect, distortion 
of the BAO rings leads to an additional effect. 

Light is deflected by large scale structure as it trav- 
els from sources in a deep survey to us, and this lensing 
leads to two sources of error in BAO. First, the position 
we assign to a given galaxy based on its redshift and an- 
gular position on the sky does not correspond to its actual 
position. When we measure the correlation function by 
assigning distance to all pairs of galaxies, we are therefore 
inevitably making an error. We place a pair of galaxies 
in a given distance bin, but they may well belong in a 
different bin. This effect tends to smooth out features in 
the correlation function so will not change the position of 
the BAO peak but will change its height and width. An 
almost identical effect acting on the CMB is well-known 
(see Ref. [2l| for a recent review) and now even included 
in the codes which compute the CMB spectrum. 

In the case of a galaxy (or QSO) survey, there is a sec- 
ond effect that does not affect the CMB. Surveys are 
generally magnitude limited, so they include only ob- 
jects brighter than some threshold. Not only does lens- 
ing change the apparent position of galaxies, but it also 
changes their magnitude p^ . [2^ : it will brighten some 
galaxies, otherwise too faint to be detected, and push 
them over threshold and de-magnify others, causing them 
to drop out of the survey. This effect, dubbed magnifica- 
tion bias, is also well-known especially in studies of the 
correlation between quasars and galaxies 0, HI] where 
it has been detected [2^, [13, H^. Its effect induces a 
non-zero contribution to the correlation function even if 
the underlying galaxy population was completely ran- 
dom [ii,!^. 

The paper is organized as follows. In Sec. Hi] we ex- 
plicitly compute the effect of the re-mapping due to the 
lensing and the smoothing it induces on the correlation 
function. In Sec. IIIII we then turn to the analysis of the 
magnification bias. In Sec. IIVI we discuss the impact of 
these two effects on BAO measurements. Finally, all the 
technicalities are gathered in the Appendices. 



II. SMOOTHING OF THE CORRELATION 
FUNCTION 

The effect we consider in this section arises because 
photons travel on geodesies and these in turn depend 
on the energy density distribution present between the 
source and the observer. This implies that the observed 
position of a source galaxy (denoted by Xo) and its actual 
position (denoted by Xs) will in general differ. The galaxy 
number overdensity at a given position x is defined by 



S3{x) 



n{x) 



(1) 



where n(x) is the galaxy number density at position x 
and h is the mean density that would be measured if all 
galaxies were uniformly distributed. 

The correlation function is defined as the two point 
function of the observed S field. The observed, lensed 
correlation function will in general differ from the un- 
lensed correlation function. What is actually measured 
is the correlation between galaxies overdensities at ob- 
served positions 

Us{r) = (,GG{r) = {5'{So)5<^{S',))\s^-S'j=r, (2) 

The unlensed correlation function on the other hand is 
the correlation between overdensities at actual source po- 
sitions ^si'f)- Since we are interested in large scales, in 
what follows we assume a linear bias model (not depend- 
ing on scales or redshift [3l|) so that galaxies trace the 
underlying dark matter fluctuations and 5^{x) = bS{x). 
This assumption in turn implies that galaxies and mat- 
ter overdensity correlation functions (resp. denoted by 
^GG and ^) are simply related by a multiplicative con- 
stant b^: in what follows we focus on the latter even 
though the results apply to the former as well. To find 
the connection between the lensed and unlensed correla- 
tion functions, we use the fact that weak lensing induces 
a re-mapping of the density field as 5{xo) = S{xs) where 
the two positions are related by the displacement vector 
field as Xq ~ Xs + Ci^s)- A Taylor expansion leads 

where we have kept terms up to second order in the dis- 
placement. It follows that (see Appendix lAl for details) 



£,ohs{r) = is{r) + 5£,sra{r), 
5^sm{r) = [(C(f)C^(f))-(C(f)C^(y))] 
'5'3 d^(r) rV^ d /$'(r) 



dr 



dr 



(4) 
(5) 



In this expression the effect of lensing appears only at 
second order because (C) = 0. 

The additional term (5^sm is due to lensing and tends 
to smooth out the correlation function. A simple way 
to understand it is to think of the BAO peak as a cir- 
cle on the sky. When the light from this circular feature 
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FIG. 1: Correction due to position re-mapping for background 
galaxies a,t z — 2.5. Solid red curve is the galaxy-galaxy 
correlation function assuming a constant bias factor 6=1. 
Green solid is the re-mapping term. Around r = 100/i^^Mpc 
it is negative while on either side it is positive. In other words 
lensing smooths out the BAO peak, as can be seen directly 
from the dotted curve which includes both contribution. 



passes through an overdense line-of-sight, it will appear 
larger so the feature will move to larger radius. Con- 
versely, it can appear smaller if the intervening matter 
is underdense. Since there is an excess of circles at the 
peak scale, most fluctuations will push a circle out of the 
peak bin and into an adjacent bin. The net result is that 
the peak will be smoothed out: the correlation function 
will be enhanced below and above the peak and reduced 
at the peak. This is encoded in equation ([5]). To further 
understand this effect, let us consider the two factors 
appearing in the lensing term separately. The zero-lag 
two-point function (C'(x)C-'(x')) is simply the variance of 
the deflection distance while {C{x)C,^ (y)) is the correla- 
tion between the deflection experienced by two light rays 
emitted a distance r from one another. If r is very small, 
then {C{^)C''{y)) ■v^ill be very close to {C{x)(^{x)), that 
is, two galaxies separated by r will have their geodesies 
deflected by the same amount. In that case, the terms in 
the first set of square brackets will cancel and there will 
be no effect: the circles on the sky will simply be shifted 
with no distortion in their sizes. However, if r is large 
enough then there will be little correlation between the 
deflections of galaxies separated by r and the second term 
will be smaller than the first. Moving on to the second 
set of square brackets, notice that the smoothing term 
can be significant where the first or the second deriva- 
tive are large. This in turn implies that the departure of 
£,oir) from ^s(r) can become non- negligible especially in 
the region of the BAO peak. 

As shown in Appendix[Bl adopting the Limber approx- 
imation all the terms appearing in equation ^ can be 




FIG. 2: Absolute value of the ratio of the correlations due to 
position re-mapping and the galaxy-galaxy correlation func- 
tion. The lensing contribution is of the order of a few percent 
in the region where the correlation function exhibits the BAO 
peak. 

recast as integrals over the dark matter power spectrum. 
The smoothing correction can then be expressed as 

Sismir) = e (y - ^2 + /a) + 7 (I - , (6) 

where 

qo2 /t4 r Ai. rXs 

r 1, z = l 

K,{k,x) = <^ Mkx), ^ = 2 (8) 

[Mkx)/ikx), * = 3 

and where W{xs,x) = 2xs(l — x/Xs) is the window func- 
tion, Xs is the comoving distance to the source, and a 
flat FRW background is assumed throughout. We inte- 
grate over the fully evolving nonlinear power spectrum 
using the procedure of Smith et al. [3^ . The results are 
depicted in Fig. [1] for simplicity in the case where all 
galaxies are at the same redshift. We see that smoothing 
occurs at the percent level, not a surprising result given 
similar conclusions in the case of the CMB. Nonetheless, 
this smoothing may need to be accounted for when an- 
alyzing future surveys. This fact is further stressed in 
Fig. [21 where the ratio between the smoothing term and 
the full correlation function is plotted for the redshift 
range that likely will be probed by future surveys. 

There are some important conclusions to be drawn 
from this analysis. First, this lensing-induced correla- 
tion is unavoidable and has nothing to do with the way 
measurements are carried out. This means that such a 
contribution will arise even if the "perfect survey" with 
an infinite limiting magnitude is carried out, simply be- 
cause geodesies are sensitive to the matter distribution 
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present between the source and the observer. Second, 
the previous derivation does not depend in any way on 
the class of objects that are being surveyed: it can be 
apphed equally well to future QSO catalogs. Third, al- 
though the corrections to the correlation function are 
larger than a percent, these corrections do not alter the 
BAO peak position. For example, at z = 2.5 this lensing- 
induced correlation weights about 1.5%, but the shift in 
the peak position is only 0.01%. BAO therefore turn 
out to be standard rulers that are very robust with re- 
spect to this effect. Finally, the dominant contributions 
to the lensing-induced terms arise from modes that are 
still in the linear regime: we carried out the integrations 
using the fully evolving nonlinear power spectrum, but 
the results do not change if the linear power spectrum 
is used. Physically this is because it is large structures 
that are most responsible for deflections; the /c-integral in 
equation ([7]) peaks at the same place the power spectrum 
does, roughly k ~ 0.02/iMpc~^, deep in the linear regime. 
An important implication of this is that the correction is 
very easy to implement. 



III. MAGNIFICATION BIAS 

We now move to consider a second, different effect that 
can also impact the determination of the correlation func- 
tion and the measurement of the BAO peak position. 
While the effect analyzed in the previous section was an 
unavoidable consequence of any metric theory of gravity 
and of the rich structure in the universe, the effect consid- 
ered here is related to the fact that actual surveys include 
only objects that are brighter than a limiting threshold. 

Weak gravitational lensing acts in two ways when such 
surveys are carried out. First, it can brighten some ob- 
jects pushing them over the threshold and de-magnify 
others causing them to drop out of the survey. Second, 
it can stretch a given patch of sky making it larger or 
smaller than it actually is, thus increasing or decreasing 
the measured number density. If no{f,z) is the num- 
ber of objects with flux larger than / per unit solid 
angle between z and z + dz in the absence of lens- 
ing, the observed number density in direction h will be 
"o[//m(", z), z]//i(n, z) where fi[{h, z), z] = - - 
7^] is the magnification expressed in terms of the con- 
vergence K and the shear amplitude 7. In the simple 
case where fi is assumed independent of z, this implies 
that n(m,n) — nQ{m)[^{h)]'^'^^~^ where s is the slope of 
the number density as a function of the magnitude m, 
s = dlnno(TO)/dm. Focusing for the time being on the 
correlation function for galaxies, it is important to keep 
in mind that the magnitude and the sign of this effect 
(over- or under-density for ^ > 1) depends on the slope 
of the luminosity function, which in turn is different for 
different classes of objects and can present redshift evo- 
lution [ill. 

Proceeding along the same lines of the analysis carried 
out by Moessner and Jain [29J on the angular correlation 



function, we find that the fluctuation in the galaxy num- 
ber density at a given position S{x) receives contributions 
from two sources: the true galaxy clustering, denoted by 
S^{x) and the variation in the number density due to 
magnification of the galaxy population induced by lens- 
ing, denoted by (5'(x). The observed overdensity is the 
sum of these two, 6(x) = S^{x) + S^x), so 

Usir) = {6ix)Sm = {S'{x)5\y)) + (5^(x),5'(y)) 

+ {S'ix)S^{y)) + {S\x)S'iy)). (9) 

where |af — ?/| — r. 

The first term in Eq. ^ is the true galaxy-galaxy cor- 
relation function. As mentioned in the previous section, 
in the linear bias approximation this term is directly pro- 
portional to the matter correlation function. 

{S<^{x)5S{y))=^GG{r)^b^^{r). (10) 

The second and third terms are the contributions to 
the observed correlation function arising from correlating 
a galaxy overdensity arising because of lensing along the 
line of sight with a true matter overdensity. The last 
term is the contribution to the correlation function due 
to two overdensities which are generated both because of 
lensing along the line of sight to the source galaxies. 

We now consider the two cross terms ((5^(x)(5'(y)). 
These arise from the correlation of overdensities due to 
clustering with an apparent overdensity due to lensing. 
To proceed further, let us recall the expression for orig- 
inally derived by Broadhurst, Taylor and Peacock jSJ 
and also used in Moessner and Jain [2^ 

6\x) ^ {5s~2)k{x), (11) 

where as before s is the logarithmic slope of the number 
counts. The convergence is determined by an integral 
along the line of sight 

where Xs is the comoving distance to the source galaxy 
and Wl{Xs,x) = xiXs - x)/Xs is the lensing geometric 
factor. Using this relation, we derive (see Appendix [C)) 
that 

c f\-/xgxi\ r^u '^ + Xsa{Xs)H{xs) j 

iGLyr) = {bH ) = Chr — ^ Igi, (13) 

where for sake of brevity we have defined the constant 
C= ^(5s-2)r!„i/2, (14) 

and the integral 

/dfc 
— P5(fc,z)Ji(fcr), (15) 

where 3\(x) is the Bessel function of the first kind. The 
factor of r in the prefactor of Eq. (|13p indicates that not 



5 



all lensing matter along the line of sight contributes to 
this cross term. Rather, only within a distance r are re- 
gions with more galaxies likely to produce a larger lensing 
signal. If there is a huge excess of matter along the line 
of sight but far from the source galaxy, the magnifica- 
tion it produces will just as likely be associated with a 
background underdensity as overdensity, so the correla- 
tion {S^S') due to it will be zero. As a result, this cross 
term will be largest for low redshift background galaxies. 
After all, the matter field is more clustered at late times, 
and the cross term is most sensitive to clustering at the 
background galaxy redshift. Fig. [3] verifies these two fea- 
tures. The cross term is indeed very small and does get 
smaller as one moves to higher redshift. 

We now turn to the lensing-lensing term in the cor- 
relation function, which is explicitly given by (see Ap- 
pendix [C| 
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where in this case 
luix) 



kdk 



Mkr)Psik,x). 



(16) 
(17) 



Here there is no prefactor of r: all structures along the 
line of sight can contribute. As a result, this term is 
largest when the background galaxies are at high red- 
shifts and can get lensed by as much structure as possi- 
ble. 

The result of the calculation of these three contribu- 
tions is shown in Fig. [3l where (2.5s — 1) = 1 is assumed 
and the integrations are carried out considering the fully 
evolving non-linear power spectrum obtained applying 
the procedure of Smith et al. [s^l ■ Four sets of curves are 
obtained, corresponding to sources at redshift of 0.5, 1, 
2, 3. From Fig. [31 we see that while the amplitude of ^gg 
and £,GL decreases for increasing redshift, the amplitude 
of C-LL increases for increasing redshift, thus contributing 
more and more to the full correlation function. 

In Fig. [5] two sets of curves are shown, obtained in- 
tegrating either over the non-linear power spectrum of 
Smith et al. [33 | or over the linear power spectrum. We 
see that non-linearities significantly contribute to the cor- 
relation function on scales that are much smaller than the 
ones where the BAO peak is located. As in the case of 
the lensing-induced correlations considered in the previ- 
ous section, this result implies that our conclusions do 
not depend on the details of the non-linear evolution of 
structure. 

It is now necessary to stress that equations (jl3p and 
(|16p respectively depend linearly and quadratically on 
the factor (2.5s — 1). This factor is crucial in determin- 
ing the weight of the magnification bias contribution to 
the correlation function of the class of objects consid- 
ered. In particular, despite the fact that the actual val- 
ues of s depend on the survey, they are smaller for galax- 
ies (s ~ [0.2,0.6]) m,!!^ than for QSO (s - [-2,1.6]) 



FIG. 3: Correlation functions at different redshifts: ^gg (red), 
£,GL (green) and ^ll (blue). In each case solid curves corre- 
sponds to galaxies at z = 0.5, dashed curves 2 = 1, dot-dashed 
curves z = 2 and dotted curves 2 = 3. For illustration pur- 
poses here we set b = 1 and 2.5s — 1 = 1 and we assume a flat 
FRW universe with = 0.3 and = 0.7. 



|2a. l37|. This means that the (2.5s — 1) factors appear- 
ing in front of the magnification bias terms are smaller 
(larger) than unity for galaxies (QSO):^ while the mag- 
nification bias is suppressed for the galaxy correlation 
function, it can actually play a significant role for the 
QSO correlation function. This in turn strengthens the 
conclusion of the previous section and points toward a 
well defined class of objects: BAO measured through 
the correlation function of galaxies represent very robust 
standard rulers, only marginally affected by lensing even 
at moderately high redshifts. By the same token, given 
the scaling of the magnification bias terms with respect 
to C, the correlation function of QSO [11] can receive 
significant contributions from the lensing-lensing term, 
especially at high redshift. If one is interested in extract- 
ing information from the lensing-lensing signal, then the 
QSO correlation function seems the way to go. 



IV. DISCUSSION 

To understand the role of lensing, recall that the mea- 
surement of the BAO bump appearing in the correlation 
function conveys information about two physical scales 
(Til [39I . The centroid of the bump is related to the sound 
horizon, that is the comoving distance that a sound wave 
can travel between the big bang and recombination. The 



^ Depending on the class of objects considered, the curves appear- 
ing in Fig. |3] and |4] need to be rescaled accordingly. 
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FIG. 4: Calculations of the (,gg, ^gl and ^ll terms carried 
out for z — 0.5 (darker curves) and z = 2.0 (lighter curves). 
The solid lines are obtained integrating over the fully evolv- 
ing non-linear power spectrum obtained using the procedure 
of Smith et al. [s^l, while the dotted ones are obtained inte- 
grating over the linear power spectrum. In the range where 
the correlation function exhibits the BAO peak non-linearities 
do not significantly contribute. Again, b = 1 and 2.5s — 1 = 1 
are assumed for illustration purposes. 



amplitude of the bump on the other hand is directly re- 
lated to the matter-radiation equality period. We have 
shown that weak lensing can affect the measurement of 
baryon acoustic oscillation through the correlation func- 
tion in two very different ways. 

First, weak lensing smooths out the acoustic peak in 
the correlation function. This smoothing is unavoidable, 
even if the "perfect survey" could eventually be carried 
out. This lensing-induced term smooths the BAO peak 



at the level of a few percent but the peak location does 
not change because of it. In other words lensing does 
not affect the measurement of the sound horizon. On 
the other hand, the measurement of the matter-radiation 
equality period (and thus of flmh^ .&]) depends on the 
peak amplitude and it can therefore be affected at the 
percent level. 

Second, weak lensing can also contribute a magnifica- 
tion bias term to the correlation function that adds to 
the unlensed term and that can potentially obscure the 
prominent feature of interest. This magnification bias 
contributes in different ways depending on the class of 
objects used to measure the correlation function. In par- 
ticular, while it can significantly affect the correlation 
function of QSO, magnification bias turns out to be only 
marginal for galaxy surveys. 

After taking both effects into account, we conclude 
that the BAO feature in the galaxy correlation function 
is a robust standard ruler. For instance, in the case in 
which all surveyed galaxies are at z — 2.5, the lensing- 
induced correlation is of order 1.5%, but the location of 
the BAO peak shifts by only 0.01%. The same cannot 
be said about possible QSO correlation function, how- 
ever, since in that case magnification bias (notably, the 
lensing-lensing term) can play a significant role. 

Finally, let us remark that the sensitivity of the results 
reported in Sec. [TTland lllll on the linear bias parameter b 
is different for the different effects. While both the lens- 
ing induced smoothing and the magnification bias would 
be affected by a scale dependence of the bias factor, the 
former would not be sensitive to a change in the over- 
all magnitude of b. On the other hand, if in the future 
evidence of "antibias" will surface (as suggested by the 
recent work of Simon et al. [Slj) this will lead to an en- 
hancement of the magnification bias. Since ^qq depends 
quadratically on b while is independent of it, hav- 
ing b ^ 0.8 would imply that the lensing-lensing term 
weights - compared to the galaxy-galaxy term - almost 
twice as much as what estimated in Sec. IIIII 
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APPENDIX A: RE-MAPPING DUE TO LENSING 

Let us compute the new term in the correlation func- 
tion due to the difference between the true position of a 
galaxy and its apparent position. Weak lensing re-maps 
the density contrast field, S, according to 



where C is the displacement field related to the deflec- 
tion angle a. For small displacement, it can be Taylor 
expanded as 



6{x + = S{x) + cs,^ + ^CCS,^, + 



(A2) 



where S^i = dS/dx'' and C are evaluated in x. It follows 
that the observed correlation function ^(r) = {6{x)S{y)) 
is given by 



^{r) + {C{x)cHy)){SA^)SAy)) 



(A3) 



e(r) + (C(x)C^(y)) 



dx^dy-i 
dx^dx^ 



e(r). (A4) 



Since = — j/*, one easily gets that d£,{r)/dx'^ 
(ri/r)^' where ^' = d^/dr and then d^£,{r)/dx^dy^ 
—d'^^{r)/dx^dx^. It follows that 



The second derivative is explicitly given by 



dmr) 
dx^dx^ 



dx^dx^ 



SijC + fifj r I — 



(A5) 



(A6) 



with = r^/r and 6ij here is the Kronecker symbol. 
Together with Eq. (jASp . this leads to the result in equa- 
tion (O. Note that |(0) = C(0), which is a consequence 
of isotropy and of the fact that {Q = 0. 



APPENDIX B: DISPLACEMENT CORRELATION 
FUNCTION 

The density field S{y) is re- mapped to the observed 
density field 6(x) with y — x + C{x) by the lensing effects. 
The displacement field C is related to the deflection angle 
a by 



(Bl) 



where Da is the comoving angular diameter distance, x is 
the radial comoving radial distance, and n the direction 
of observation [so that choosing the coordinate system 
orientation with the z axis aligned along the line of sight 
X ~ (£)^(x)n, x)]- This means that we assume that the 
effect of gravitational lensing is to shift the position of 
the galaxies perpendicularly to the line of sight and that 
we are neglecting the (small) effect on the redshift and 
thus on the radial coordinate. 

The deflection angle is related to the gravitational po- 
tential $ integrated along the line of sight by 



d{n,x) 



~S{x) = S{y)=S{x;+0, 



(Al) 



DaIx-x') 

Da{x) 



(B2) 
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where V_l is the gradient in the 2-dimensional plane per- 
pendicular to the line of sight [13, |4l| . It follows that 

dh, x) = W{x, x')^±'^[DA{x')n, X'W, (B3) 
with the window function 

Wix,x') = 2DA{x~x'), (B4) 

which reduces for a universe with Euclidean spatial sec- 
tions to W{x, x') = 2x(l " x'/x)- 

The gravitational potential is expanded in harmonic 
space with the Fourier convention 

^^-p^<^>(fc,77)e''=■^ (B5) 



and we define the power spectrum 



($(fc,77)$*(fc',r;')) = P*(fc;?7,r;')(5('Hfc-fc'), (B6) 



where rj is the conformal time so that the unperturbed 
geodesic equation is x = Vo ^ V- 

Using Eq. (|B3[) and the Fourier decomposition, the dis- 
placement correlation function is given by 



{CmH.y)) = w{x.,xi)dxi W{xy,x2)dx2 J^^^^,), ^^^^3% VYiVi,($(fci,xi)$(fc2,x2))e'('-^-^-"--^), 



where X and Y are the equation of the null geodesies re- 
lating the galaxies respectively in x and y to the observer. 
We do not state here whether x — y or not. In spher- 
ical coordinates of the position x (resp. y) is ixxiftx), 
fix being a spacelike unit vector and X (resp. Y) given 
{Xi,nx) with xi=T]o- V- 
We use Eq. (jBSp to express the field correlator and 

integrate over k2 to get 

iCm'm - rW{xx,Xi)dxi rV(x„X2)dx2 
Jo Jo 

/d'^fc 
^^-^fclfc^P<i.(fc;xi,X2) 

where we have decomposed the wave-vector as k — 
{k±, fc||). On small scales (9 <C 1), we can make the Lim- 
ber approximation (23l . l4l| which exploits the fact that 
the main contribution to the signal comes from trans- 
verse modes so that P<s>{k;x,x') — P<s>{k±;x,x')- This 
allows us to integrate over fc|| to get 27r^'^^^(xi — X2) and 
then on X2 to get 

fXM 

iCmHy)) = / W{xx,x')W{xy,x')dx' (B9) 
/ ^klkiP^{k^;x')e''^-'^\ 

with i?x = Xj_ — Yj_ and where xm = niax[xa;, Xj/]- 
order to evaluate Eq. (IA6|) . we need to compute both 



{C{^)Ci{y)) and {'fi('-{x)fj(^{y)), paying attention to the 
fact that X may or may not be equal to y. 

Let us first concentrate on {(^{x)(i{y)). We choose the 
orientation in the plane perpendicular to the line of sight 
such that k± ■ R± = k±xcos(p. If x = y, the integra- 
tion over just gives 2tt. \i x ^ y, we have to inte- 
grate exp(ifc_LX cos 1^9), which can be easily seen to give 
27rJo(fc_Lx) by using the series expansion of the complex 
exponential in terms of Bessel functions, 

00 

g^xcosv = Jo(a;) + 2^^"cos(n^)J„(x). (BIO) 

n=l 

We conclude that 

1 rXM 

(C(x)C,(y)) = — / M^(x.,x)W^(Xy,x)dx (Bll) 

/{io(fcx)}''*^'^'^'^^'^'^'^'^°'^ 

Let us now turn to the second term {Q{x)Cj> {y))fif.y 
Again, we choose the orientation in the plane perpendic- 
ular to the line of sight such that k±-R± — fcj^x cos tp. We 
have term fc^ is now replaced by {k]_ ■ fi)^ = fc^ cos^ (p. 
If X = y, the integration over (p of cos^ (p just gives tt. 
If X ^ y, we have to integrate cos^ iy9exp(ifc_LXCos(p). 
Again, we use the expansion (jBlOp and only two terms 
contribute to the integration over ip: Jo{kx) cos^ ip gives 
7rJo(fcx) and for n = 2, 2P cos2(pcos^ p gives — 7rJ2(fcx)- 
We conclude that 
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1 

47r 



Mkx)~J2{kx) 



P^ik;x,x)k^dk, for ;! | 



(B12) 



It follows that Eq. (jA6[) reduces to 

i{r)^ar)+Ak' + Br(^) , (B13) 
r \ r J 



with 



1 r r^M 

A = — k^dk 

27r J .,0 



dxP'i'ik;x,x)W{xx,x) 



X [W{x.,x)~W{xy,x)Mkx)], 



(B14) 



^ = ^^'^^ J dxP^{k;x,x)W{x.,x) 

X [W^(X.,X) - W{xy,x){Mkx) - J2(fcx))lB15) 

To finish, one should give an expression for P^. The Pois- 
son equation on sub-Hubble scales, A<i> = inGpa^S = 
3r2„ioi?oi5/a implies that 



P'i>{k;x,x') = Ps{k;v,v') 



V 2fc2 



.9(0)^ 



where g ^ D/a with D being the growth factor of 
the density perturbations, 5{k,ri) = D(r])6init{k), and 
Ps{k;0) = Ps{k) the density contrast power spectrum 
at redshift 0. 



APPENDIX C: GALAXY-LENSING AND 
LENSING-LENSING CORRELATIONS 

In this appendix, we give explicit expressions for the 
correlators needed for the computation of the magnifica- 
tion bias. 



1. Galaxy-lensing term: {Srv^Sri^) 

First, let us consider the galaxy-lensing cross correla- 
tion function. It corresponds to the second [and third] 
term in equation Going back to the second term of 
the sum we have 

{6nl{xo)5n'*{yo)) = (5s ~ 2)(<5n?(xo)«;* (yo)), (CI) 

where 5n is the fluctuation of the galaxy number density 
and 5n} has been related to the convergence n by Eq. pTjl . 
The l.h.s correlator is given by 



{5n{{xli)Hi*{yo)) 



XLivo) w 

rflraH^b I dx—{S{xh)S*{y{x))) 



_ ^fYi.H(\ b 



XL(yo) w 

dx— C[r(x)], (C2) 



where f{x) — \xi) — y(x)\ so that r(x) the distance be- 
tween the source at afg and the point y(x) which is run- 
ning along the line of sight between the observer and the 
source point at j/q. Wl is the lensing geometric factor. 



Wl{xl,x) = xiXL ~ x)/XL, 



(C3) 



and £_{r) is the unlcnscd (or "true") matter correlation 
function 



C[Kx)] = 



k^dk 



3o[kr{x)]P{k)^^{x). 



(C4) 



Equation (|C2p tells us that the correlation between a 
true matter overdensity and a lensed one is given by an 
integral over the line of sight of the correlation function 
relative to the true overdensity and the point running 
along the line of sight weighted by the usual lensing 
factor. We can actually work on the above expression 
much further in the spirit of Limber's approximation. 
Letting for brevity (SncSriL) = {Sn^{ocQ)Sn'^ (yo)) and 
C = (5s — 2)3rimiJo/2 and remembering the dependence 
on the sources' positions of the two terms we have 

(<5nG<5ni) = Ch \ dx— -^e^''^^^^ Ps{k, z) 
Jo a J {2ttY 



(27r)3 



a 



X e 



i[kiXGSi+k2XGS2 + k:-i{xG~X)] 



(C5) 



where in the last step we have picked the coordinate sys- 
tem so that ?7o — (0, 0, xl), = Xg(^i, 6*2,1) and the k 
coordinate system so that fcs is aligned along j/q.^ 

So far we have been considering the completely general 
case in which xg 7^ Xl- For simplicity let us now set 
XL = XG = Xs- Then 



fXa(!7o) 

{SncSnL) ^ Cb I dx 



X e 



J (27r) 

i[kiXsdi + k2XsS2 + k3{x 



^^-^Psik, z)— ^ 

x)] 



a(x) 



(C6) 



^ The result will not change if instead we pick yo = Xl{ 
and xo = xg(0.0,1). 



L,e2,l) 
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Using the fact that 

WLexp[ik3{xs - x)] = -i{xlXs)dk3 exp[iA:3(Xs - x)], 

(C7) 

and then integrating by parts over leads to 



where we have used the fact that for any function 
of the modulus only dk-^fik) = (fc3/fc)d/(fc)/dfc = 
{k3/k±)dk^f{k) with fc^ = k'j_+kl and where, again, we 
have used that k^ exp[ifc3(xs — x)] = '^'^x 6xp[i/c3(xs — x)]- 
Now in the Limber's limit in which the k^ dependence of 
the power spectrum is neglected, we can integrate over 
^3 to get 27r(5(^) (xs — x) (where 5^^^ is the derivative of 
the Dirac distribution). Integrating over x then leads to 



{dncSnL) 



Cb 



d^k 



{l+ClsXsHs) 



QsXs J (Stt) 

X 9fc,P5(fci,Xs)e'['^^^'^+'^^^'^i, (C9) 
I 



with Os = a(xs) and where we have used the fact that 
da/dxs = —O'lHs with = H{as)- We integrate over 
the angle with the use of Eq. (jBlOP to get (27r) Jo(r) with 
r — XsO and then make an integration by part to get 



{SncSriL) — Cb 



1 



ClsXs 



dk 
2^ 



Psik,z)Ji{kr). 



(CIO) 



This leads to the result in Eq. ([T 

2. Lensing-lensing term: {Sn{Sn2) 

The lensing-lensing term 

{6n[{x'o)6n'*{yo)) = {5s ~ 2f {k{xo)k* (yo)) ,{011) 

is directly related to the convergence correlation function 
and can be computed in a similar way as the galaxy- 
lensing term. We start from the expression 



{K{xo)K*{yo) 



rXiiSo) r 

/ '^^^ / 

Jo ^0 



Xi(xi(^o) 

«X2 



Xi) X2 (X2(27o) - X2) 



Xi(a^o)a(xi) X2(j/o)a(X2) 



{S{x{xi))S{y{x2)))- (C12) 



r 



Indeed, when the two sources lie on the same line of sight, 
it reduces to the variance of the convergence. However, 
in general the two lines of sight - parameterized by a^(xi) 
and y{x2) ~ are not identical. In this case, defining r' = 
\x{xi) — y{X2)\, it is possible to recast Eq. (|C12p in terms 
of the correlation function £,{r'). 

We use the same conventions and notations as in the 
previous paragraph. We choose the coordinate system 
such that y(x2) = (0,0, X2) and x{xi) = Xi(^'i,^2,l) 
and we orient the fc-space coordinate system so that ^3 
is parallel to y{x2)- It follows that Eq. (|C12p takes the 
form 

Jo Jo a(Xi) a(X2) 

X [ -^P5(fc,z)e*['=ixi«i+'=«i«^+'=^(xi-X2)]^Ci3) 
J (^tt) 

Again, we use the Limber's approximation so that 
Ps{k,z) ~ Ps{k±,z). In this limit, the integral over ^3 



simply gives (27r)(5(^^ (xi — X2) and 



{SnLSriL) ^ dx 



Wl{x) fd'k, 

^ a'{: 
ix(kiei+k2e2) 



■Ps{k,x) 



lo aHx) J {-^^f 
X e'^^'"-'"-' ' ^ (C14) 

To finish, we integrate over the angle in the /c^ -plane to 
get 



{5nL5nL) = CM dx 
Jo 



Wl{x) 



n 2 



a(x) 



Ill{x), (C15) 



with 



Iu{x) 



kdk 



Jo{kr)Ps{k,x), 



(C16) 



which is exactly the result given in Eq. p^ . 



